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ABSTRACT 

We study the transition of nuclear matter to strange quark matter (SQM) inside neutron stars (NSs). 
It is shown that the influence of the magnetic field expected to be present in NS interiors has a dramatic 
effect on the propagation of a laminar deflagration (widely studied so far), generating a strong acceleration 
of the flame in the polar direction. This results in a strong asymmetry in the geometry of the just formed 
core of hot SQM which resembles a cylinder orientated in the direction of the magnetic poles of the NS. 
This geometrical asymmetry gives rise to a bipolar emission of the thermal neutrino- antineutrino pairs 
produced in the process of SQM formation. The vv annihilate into e + e~ pairs just above the polar caps 
of the NS giving rise to a relativistic fireball, thus providing a suitable form of energy transport and 
conversion to 7-emission that may be associated to short gamma ray bursts (GRBs). 

Subject headings: Stars: Neutron - Gamma Rays: Bursts - Instabilities, Magnetohydrodynamics: MHD 
- Stars: Magnetic Fields, 



1. INTRODUCTION 

The transition to SQM is expected to occur in NSs if 
SQM has lower energy per baryon than ordinary nuclear 
matter Bodmer 1971, Terazawa 1979, Chin and Kerman 
1979, Witten 1984. This has strong consequences in the 
astrophysics of compact stars since, if SQM is the true 
ground state of strongly interacting matter rather than 
56 Fe, compact objects could be strange stars (SS) instead 
of NSs. Some tentative strange star candidates are the 
compact objects associated with the X-ray bursters GRO 
J1744-28 Cheng et al. 1998, SAX J1808.4-3658 Li et al. 
1999a, and the X-ray pulsar Her X-l Dey et al. 1998. In 
addition, the observed high and low frequency quasi peri- 
odic oscillations in the atoll source 4U 1728-34 have been 
shown to be more consistent with a strange star nature 
Li et al. 1999b. A number of different mechanisms has 
been proposed for NM to SQM conversion Haensel and 
Madsen eds. 1991 inside the star. All them are based on 
the formation of a "seed" of SQM inside the NS. A pos- 
sible mechanism is the so-called strangelet contamination, 
where a seed of SQM from the ISM enters a NS and con- 
verts it to a SS. In another possible scenario, a seed of 
SQM forms in the core of a NS as a result of the increase 
of the central density above the critical density for decon- 
finement phase transition. A way to do this is through 
mass accretion onto the NS in a binary stellar system. In 
a third mechanism, a seed of SQM may naturally form in- 
side a newly born neutron star from a core-collapse super- 
nova explosion after a deleptonization timescale Lugones 
and Benvenuto 1998, Benvenuto and Lugones 1999. No 
matter which of these mechanisms are actually triggering 
the NS to SS transition, once the first seed of SQM is pro- 
duced inside the NS, it will propagate as a combustion 
swallowing neutrons, protons and hyperons. The transi- 
tion to SQM inside the burning front has actually two 
stages. Deconfinement driven by strong interactions first 
liberates quarks confined inside hadrons. The just decon- 



fincd quark phase has a certain finite strangeness due to 
the presence of strange hadrons in NS matter. However, 
the composition of the just deconfined phase (with u,d 
and s quarks) is not in beta equilibrium and consequently 
chemical equilibrium is reached by weak interactions. It is 
in this phase that a large amount of energy is released and 
then v's are copiously produced. 

The type of diffusion-driven combustions studied so far 
by several authors Olinto 1991, Heiselberg et al. 1991 ac- 
tually correspond to laminar deflagrations (slow combus- 
tions). Whether the conversion process remains forever 
as a deflagration, either laminar or turbulent; or jumps 
to the detonation regime, thus driving an explosive tran- 
sient Horvath and Benvenuto 1988, Benvenuto and Hor- 
vath 1989, Lugones et al. 1994 has been debated in the 
literature. The situation is closely analogue to the much 
more studied thermonuclear combustions leading to type 
la supernovae. In addition to the lack of detailed studies 
on the mode(s) of combustion (laminar/turbulent defla- 
gration or detonation), there is (to the best of our knowl- 
edge) no calculation of the influence of ubiquitous mag- 
netic fields expected to be promptly generated or already 
present in the NS as a fossil. We shall consider hereafter 
an initially laminar deflagration in presence of a £?-field. 
As stated, previous works Olinto 1991, Heiselberg et al. 
1991 have calculated the velocity of the laminar deflagra- 
tion by considering the diffusion of s-quarks as the main 
agent for the progress of the conversion. Given that the 
combustion is idealized to happen near the T — limit 
(which is small when compared with the chemical poten- 
tial of quarks), these diffusion-limited "cold" models are 
reasonable for this purpose. The general result, which also 
holds when temperature corrections and full non-linearity 
are considered, is that the laminar velocity of the front 
is relatively slow (vi am < 10 4 cm/s) although this speed 
may be uncertain by several orders of magnitude (Olinto 
1991). This velocity is a direct consequence of both the 
timescale for weak decays that create s-quarks (10 _8 s) and 
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the physics of the diffusion process. 

2. INSTABILITIES AND ASYMMETRIES 

A nuclear flame that starts as a laminar deflagration 
propagating outward rapidly enters the wrinkled flamclct 
regime due to the action of several hydrodynamic insta- 
bilities, such as the Landau-Darrieus (LD) and Rayleigh- 
Taylor (RT). For the largest scales (those ^> than the 
thickness of the flame) , the RT instability dominates over 
LD (see Ghezzi et al. 200f for details). It is possible to 
identify distinct regimes in the deflagration stage, and we 
will apply the fractal model of combustion Timmes and 
Woosley 1992, Niemeyer and Hillebrandt 1995, Niemeyer 
and Woosley 1997, Ghezzi et al. 2001 for all those regimes. 
The wrinkled surface A behaves like a fractal with A oc 
R D , where D is the fractal dimension of the surface, with 
2 < D < 3, and R is the mean radius of the wrinkled sur- 
face Filyand et al. 1994. Numerical simulations Filyand 
et al. 1994 and laboratory experiments involving different 
gas mixtures Gostintev et al. 1988 show that the fractal 
growth actually increases the velocity of the combustion 
front because of the change in the transport mechanism 
from a laminar to a fully turbulent burning. The burn- 
ing of nuclear matter to SQM has been studied so far in 
the absence of any £?-field. However, there is an ample 
consensus that strong magnetic fields should be present 
at the center, as indicated by the observations of surface 
pulsar fields ( B > 10 n ~ 12 G ). Even though there is no 
detailed information on the fields of non-pulsating NSs, it 
seems reasonable to assume a field at least as strong as the 
one at the surface of radio pulsars inside the object. 

The fractal model presented in Ghezzi et al. 2001 for 
the thermonuclear burning of a white dwarf is applica- 
ble to the present case since all the basic ingredients and 
the key points of the physics are similar. For the sake 
of simplicity, we shall assume hereafter a dipolar geom- 
etry and thus the flame propagating in two particularly 
representative directions, one parallel to the B-field lines 
(in the polar direction) and the other one perpendicular 
to them (in the equatorial direction). Within the fractal 
description, the effective velocity of the flame is given by 
v f — vi am (L/l) D ~ 2 where L and I are the maximum an 
minimum length scales of perturbations unstable to RT . 
The characteristic velocity of the RT growing modes must 
be > vi ami i. e. Innxil) — vi am , where urt is the inverse 
of the character istic RT t ime. In the absence of a B-ficld 
n RT = {l/2n)y/gkAp/2p and so I = 4npvf am /gAp. How- 
ever, in the equatorial direction, the presence of B is essen- 
tial in modifying the dispersion relation for the RT insta- 
bility which reads urtb = l/2Tr v / gk(Ap/2p - kB 2 /4ngp) 
where k — 2n/\, A is the wavelength of the perturba- 
tion, and Ap = p u — pd is the density difference be- 
tween the upstream and downstream parts of the flame 
front Ghezzi et al. 2001. Taking into account the ve- 
locity of the flame, the minimum cut-off length is given 
by the condition l e riRTB(le) — vi ami from which we de- 
rive the minimum scale in the equatorial direction l e = 
8tt(B 2 /8tt + pvf am /2)/gAp. So, from the above definition 
of the fractal velocity {L/l e ) D - 2 and 

1 We should note that, although the present work is based solely 
magnetic RT instability (Jun et al. 1995) reinforces the results of as; 
the non- linear regime. 



v p = vi am (L/l p ) D ~ 2 . Note that L is not modified by B 
and has the same value in both directions. Then, the ratio 
between the equatorial and polar velocities is Ghezzi et al. 
2001: 

^v p /v e = [l + B 2 /(4npv 2 am )] D - 2 . (1) 

That is, although the magnetic pressure B 2 /(8w) is not 
relevant, the B-ficld quenches the growth of RT instabil- 
ities in the equatorial direction acting as a surface ten- 
sion, while it is innocuous in the polar one where in av- 
erage we have v p x B = . More specifically B modifies 
the minimum RT instability scale and since the turbu- 
lent flame velocity is related to RT-growth this results in 
a different velocity of propagation along each direction. 
Since the density of the products of combustion (SQM) is 
comparable to the fuel (nuclear matter) p ~ 10 15 <?/cm 3 ; 
and using the values for vi am found in previous works 
Olinto 1991, Heiselberg et al. 1991 (~ 10 4 cm/s), we 
find that B 2 /(4tt) > pvf am for relatively low values of 
B(i. e. ~ 10 12 ~ 13 G), so that for D w 2.5 (which is in good 
agreement with numerical studies, Blinnikov et al. 1995), 
£ scales linearly with the field to a good approximation. 
So, 

£ ~ 10 x (B/10 13 G)(10 15 gcm- 3 /p)^ 2 (10 i cm S - 1 /v lam ) 

(2) 

and large asymmetries can be produced even for moderate 
values of B 1 . Next, we evaluate the time that the polar 
front needs to reach the surface (t p = R/v p = R/(^v e )): 

t p ~ 10s x (R/10km)(10 13 G/B)(p/10 15 gcm- 3 )^ 2 . (3) 

Note that eq. (3) gives actually an upper bound to r p be- 
cause we have approximated v e by vi am . An asymmetry in 
the ^-emission will be possible only if r p < r<j, that is, if r p 
is smaller than the typical diffusion timescale Td through 
nuclear matter in the equatorial direction (jd is at least 
30s, Pons et al. 1999). Only in this case we can assume 
that the lateral sides of the cylinder are opaque and almost 
all the vv are emitted through the polar caps. Otherwise, 
we will still have a strong asymmetry in the shape of the 
SQM region but this will not yield an asymmetry in the 
^-emission structure. The minimum B for which t v < Td 
is satisfied is obtained from Eq. (3) 

Bmin = 3 x 10 12 G{p/l0 15 gcm- 3 ) 1/2 {R/10km). (4) 

In other words, for B < B min we recover the isotropy of 
the neutrino emission (see below). 

When the polar front reaches the stellar surface, the 
equatorial front will have travelled only 5 = R/£ = 1km x 
(R/10km)(10 13 G/B)(p/10 15 gcm- 3 ) 1 / 2 {vi am /10 4 cms- 1 ), 
and therefore the SQM region will begin to emit its neu- 
trino content mainly through a cylinder aligned with the 
poles with a radius <5 ~ 1 km and a height 2R ~ 20 km (see 
Fig 1). Notice that if B 10 13 G, 8 is very small and the 
total energy released by the SQM conversion is negligible 
in this asymmetric stage. An estimate for the ratio of the 
neutrino flux in the equatorial and the polar directions 
is found by comparing the flux through the opaque NS 

a linear analysis, 3D numerical simulations of the development of the 
netry here reported and also reveals a tendency for its amplification in 
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matter in the equatorial direction with the free streaming 
flux of neutrinos emitted from the polar caps, namely 

F e /FP = [(-cA„/3) dejdr] / [«„] « A„/(3Ar) (5) 

where Ar is the difference between the radius of the 
star and the equatorial radius of the SQM core. \ n = 
3 x 10 3 cm(p nuc / p){10MeV/ E v ) 2 is the mean free path of 
neutrinos diffusing from the surface of the SQM core in the 
equatorial direction Shapiro and Teukolsky 1983. This A„ 
is mainly due to neutral-current elastic scattering off neu- 
trons. For typical values of A„ and a highly asymmetric 
SQM region' we have F e /F p ~ 1CT 2 - 10" 4 . Thus, once 
the polar front reaches the surface, almost all u's are re- 
leased through the polar caps. This is due to the fact that 
the flux in the equatorial direction is quenched due to the 
high A„ while vv 's in the polar surface escape freely. 

The neutrino transport inside the SQM cylinder can 
be described by: detot/dt = d / dx[{\sQMc/2>) de v /dx], 
with the boundary conditions de v {Q,t)/dx = He u (0,t), 
de u (l,t)/dx — He u {l,t) and the initial condition 
e v (x, 0) = e„ , being e v (x,t) the neutrino energy density, 
etot{x,t) the total energy density of SQM, I — 2R the 
height of the cylinder, and H = 3/Xsqm- For the sake of 
simplicity, we assume that neutrino transport is driven by 
electron-neutrino pairs, so that Xsqm = 1-3 x 10 3 cm/T 1 2 1 
Haensel et al. 1991. Although muon and tau neutrinos are 
in fact important, this simplification does not essentially 
modify the picture and allows a simple solution for the 
problem. The total internal energy released in the con- 
version of a whole NS into SQM is £/ ~ few x 10 53 ergs 
Bombaci and Datta 2000. So, a rough estimate of the ini- 
tial temperature of the just formed SQM will be T\\ ~ 3 
(in units of lO 11 ^). From the integration of the equations 
above we obtain the temperature in the polar cap surface 
as a function of time T p u(t), and the corresponding neu- 
trino luminosity 

L v (t) = 1.2 x 10 53 erg/s [T pll (t)/3] 4 (5/lkm) 2 . (6) 

The emitted vD's annihilate into e+e~ pairs with rather 
low efficiencies. Following the procedure of Haensel ct al. 
1991 we have calculated the luminosity in e + e~ for this 
cylindrical configuration (see Haensel et al. 1991 for the 
result in spherical symmetry) 

L pair (t) = 2.1xl0 52 erg/s [T pll (t)/3] 9 {S / lkm) 2 {R/ 10km). 

(7) 

These luminosities are shown in Fig. 2. We find that 90% 
of the e + e~ pairs are injected inside small cylinders lo- 
cated just above the polar caps (with radius 6 and height 
0.4R) in a timescale of r* ~ 0.2s almost independently on 
the initial temperature. 

3. APPLICATION TO GRBS 

GRBs appear to fall into at least two distinct categories, 
namely the short duration bursts (~ 0.2 s) and the long 
duration ones (« 20 s) (for a review see Piran 2000). Pre- 
vious works have explored the idea that the conversion of 
nuclear matter into SQM in NSs may be an energy source 
for GRBs Alcock et al. 1986, Ma and Xic 1986, Haensel 
et al. 1991, Cheng and Dai 1996, Bombaci and Datta 

2 We should note that according to recent calculations of proto-NS 



2000, Wang 2000, Ouycd 2002. These models addressed 
spherically symmetric conversions of the whole NS giving 
isotropic 7- emission. We show here that if a conversion to 
SQM actually begins near the center of a NS, the presence 
of a moderate magnetic field B (~ 10 13 G) will originate a 
prompt asymmetric 7-cmission which will be observed as 
a short, beamed GRB. Given that all pairs contributing to 
Eq.(7) will annihilate into 7's, we are lead to the conclu- 
sion that an asymmetric short GRB can be generated by 
this process. The energy injected in 7-rays is E 1 > 10 51 
erg, comparable to the ones derived for long GRBs. The 
asymmetry crucial for this conclusion is controlled by the 
quotient of the field B to the laminar velocity vi am (Eq.l). 
Had the latter been larger by a few orders of magnitude no 
significant asymmetry would be produced by realistic mag- 
netic fields. It should be remarked that we have not shown 
that the turbulent deflagration will be the actual propaga- 
tion mode, but rather assumed its realization all the way 
down to the stellar surface. Actually, a possible transi- 
tion from a deflagration to a detonation mode (DDT) (in 
close analogy with the DDT in white dwarf combustions), 
may be achieved Horvath and Benvenuto 1988, Benvenuto 
and Horvath 1989, Lugones et al. 1994 depending on the 
behavior of the microphysics. However, if a substantial 
asymmetry has already been produced when the DDT is 
achieved, the polar flow will detonate before the equatorial 
one. 

Basic requirements for associating asymmetric SQM 
burnings to (short) GRBs arc the absence of a large 
amount of matter above the polar regions and a rate per 
galaxy high enough (~ 10 _3 yr _1 ) so as to compensate 
the collimation of the emission. It is unclear whether 
both can be achieved, but in any case it is clear that fur- 
ther studies of the hydrodynamics and statistics of the 
events are desirable. However, and quite independently 
of the specific astrophysical setting, the result of NM to 
SQM conversion in presence of B ~ 10 13 G will produce 
an asymmetry in the neutrino emission and an enhanced 
annihilation to 7's. However, the observational outcome 
and the rate of the events will depend on the type of sys- 
tem in which the NS is being burned. For example, in 
transitions produced by accretion induced collapse of a 
white dwarf, the energy is deposited in a region that con- 
tains a low baryon loading. So, this collimated emission 
has an advantage over the isotropic one, where B = 0, in 
that the total baryon loading "seen" by the beamed burst 
(~ S 2 /(4R 2 ) = Trpv 2 am /B 2 ~ 10- 2 - 10- 3 of the total) 
is low enough to allow for a relativistic expansion of the 
fireball with very high Lorentz factors and is sufficient to 
explain apparent burst luminosities L 7 up to more than 
~ 10 52 ergs/s for burst durations of t 1 « 0.2s. NSs in 
Low Mass X-ray Binaries (LMXBs) are also likely candi- 
dates for the conversion to SQM because mass accretion 
may raise the central density of the NS above the critical 
density for transition. These NSs have "weak" magnetic 
fields in their surface (B = 10 8 ~ 9 G) which are generally 
understood as the result of mass accretion. Since this af- 
fects in principle only the field in the surface, the field 
in the interior of the NS would be in the range of values 
that leads to an asymmetric short GRB. Note also that 
if the emission is beamed it would be difficult to explain 

evolution, convective instabilities occur within the first second of NS 
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the observed rate of GRBs by conversions in binary sys- 
tems alone. For transitions in newly born proto-NSs 2 , the 
presence of the supernova envelope makes necessary to ex- 
plain how the 7-emission produced by the central proto-NS 
traverses the young expanding ejecta (one clue may have 
been provided by GRB980425/SN1998bw association). In 
any case, we may expect emission of quasi-thermal X-ray 
flashes with rates per galaxy r/lO^yr' 1 , being r\ the un- 



known fraction of proto-NS with fields B ~ 10 G from 
the events. Other consequences for the ejected envelopes 
may include asymmetries of the envelope itself due to the 
non-isotropic injection of energy. We plan to discuss these 
issues thoroughly in a future work. 
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Fig. 1. — A moderate field (B ~ 10 13 G) in the interior of a NS being burned to SQM is able to generate an asymmetric velocity field in the 
phase change front which, depending on the exact value of B, may reach values 3> 1. So, a prolate SQM fluid is the actual geometry achieved 
by the conversion in presence of such non-zero B. This affects the leaking of thermal neutrinos produced in the transition, funnelling them 
into opposite polar cones. 



Fig. 2. — Neutrino and pair luminosities for B = 10 13 G. The timescale on which 90 % of the pairs are injected is 0.2s. 
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